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Abstract. In a single molecule assay, the motion of a molecular motor is often 
inferred from measuring the stochastic trajectory of a large probe particle attached 
to it. We discuss a simple model for this generic set-up taking into account explicitly 
the elastic coupling between probe and motor. The combined dynamics consists of 
discrete steps of the motor and continuous Brownian motion of the probe. Motivated 
by recent experiments on the Fi-ATPase, we investigate three types of efficiencies both 
in simulations and a Gaussian approximation. Overall, we obtain good quantitative 
agreement with the experimental data. In particular, we clarify the conditions under 
which one of these efficiencies becomes larger than one. 
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1. Introduction 

Molecular motors are protein complexes of the size of nanometers that convert chemical 
energy into mechanical motion [U [2]. Operating in an aqueous solution they exhibit 
stochastic dynamics and energetics due to the influence of thermal fluctuations. 
Unbalanced concentrations of the molecules providing chemical energy as input cause 
the motor proteins to operate under nonequilibrium conditions which induces a rectifled 
motion with non-zero average velocity. Consequently, molecular motors are often 
modelled using Langevin, Fokker-Planck or master equations. The so called ratchet 
models combine continuous diflFusive spatial motion with stochastic switching between 
diflFerent potentials corresponding to diflFerent chemical states j3[ Sj. Alternatively, 
transitions among a discrete state space governed by master equations provide another 
possibility to model molecular motors [SJ [6l [H El [9] . 

A quantity of general interest is the eflBciency of such stochastic machines (iQl [HI 
[T2l [T3l [T4] . For motor proteins, different kinds of efliciencies can be deflned depending 
on whether one focuses on the work against an external force, i.e., thermodynamic 
eflBciency or whether work against viscous friction is also taken into account like in the 
Stokes or generahzed efliciency [ISl [lEl [HI [H [IHl [20] . 

Experimentally, the properties of motor proteins can be investigated in single 
molecule experiments by attaching probe particles of the size of micrometers to the 
motor protein and by observing the trajectories of the probes. Additionally, such probes 
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allow to exert forces on these motor proteins [22[ |23l 12101 . Literally speaking, in these 
assays one cannot observe the motion of the motor directly but rather has to infer its 
properties from analyzing the trajectory of the probe particle. Generically, some elastic 
linker couples these two elements. Inferring properties of the motor protein requires to 
consider the interaction effects that depend on the linkage between motor protein and 

probe [ZSlEelEzlEslEniEQlE^ 

In the present paper, we discuss a minimal hybrid model for such a motor protein 
assay that includes this elastic link explicitly. The motor protein and the probe will be 
modelled as two degrees of freedom moving along a spatial coordinate. In particular, 
we investigate different kinds of efficiencies used previously to describe the energetics of 
molecular motors and compare our results quantitatively to recent experiments of the 
rotary motor protein Fi-ATPase [33l EH [35l 06]. Previous theoretical modelling of Fi^ 
ATPase using a discrete state model as well as a ratchet model assuming the probe to 
stick directly at the motor has especially focussed on the dependence of the rotational 
behaviour on friction, external forces, nucleotide concentrations and temperature as 
well as on chemical and thermodynamic efficiency and the fluctuation theorem f37l 138] . 
Detailed modelling of the rotary mechanism and the involved subunits can be found in 
[391 So]. 

This particularly well studied molecular motor consists of three a and three ^ 
subunits arranged around a central 7 shaft |41j- Binding and hydrolysis of an ATP 
molecule at a /3 subunit drives a rotation of the 7 shaft of 120° [2T] which has been 
observed to consist of two substeps of 90° and 30° |42]. An external torque exerted on 
the 7 shaft (as experimentally done in [23] or by the Fq part within the cell) induces 
ATP synthesis. Coupled to the membrane embedded Fq part, Fi^ATPase provides 
ATP for further hydrolysis reactions therefore being an important part in the energy 
transfer of cells. Experimental observations of the Fi-ATPase in the hydrolysis direction 
include the measurements of different kinds of efficiencies. The Stokes efficiency, a Stokes 
efficiency confined to single jumping events and the thermodynamic efficiency, especially 
at stall conditions, have been investigated [43l SU |34] . These experiments led to values 
for the Stokes efficiency and the thermodynamic efficiency of almost 1 suggesting that 
the Fi-ATPase can use almost the complete chemical energy either to drive the probe 
through a viscous medium or to perform work against an external force. Recently, a 
measure of the efficiency that takes explicitely care of fluctuations was introduced [33j . 
The definition of efficiency used there also provided values close to 1 for the examined 
parameters. Our analysis will show that the latter efficiency can easily reach values 
larger than 1. 
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Figure 1. Schematic representation of the motor protein (blue) with attached probe 
(red). The instantaneous distance between motor protein and probe is denoted by 
y(t). The probe moves along a continuous spatial coordinate x{t) an is subject to an 
external force /ex- With transition rates k^{nj^x) the motor protein jumps at times 
Tj between discrete states nj separated a distance d. The load sharing factors ^+ and 
0- indicate the position of an underlying unresolved potential barrier relative to the 
potential minimum. 



2. Hybrid-Model 

2.1. Single molecule dynamics 

The one-dimensional model we will use to describe a molecular motor with an attached 
probe particle consists of two degrees of freedom representing the motor protein at 
position n{t) and the probe at position respectively, see figure [H For a rotary 

motor like the Fi-ATPase, the rotary motion is mapped to a linear one for simplicity. 
Both constituents are linked via a harmonic potential 

V{n^ x) —{n — xY (1) 

with spring constant where we have included a possible rest length of the linker into 
the definition of x. 

The motion of the probe particle is described by an overdamped Langevin equation 
with friction coefficient 7 and constant external force /ex, 

i^W = (-^-/ex)/7 + CW, (2) 

including the random force 7C(^) that the solvent exerts on the probe. The thermal 
fluctuations are assumed to be Gaussian white noise with zero mean and correlations 
{C{'ti)C{h)) = 2/cBT5(ti— ^2)77 where /cb is Boltzmann's constant and T the temperature 
of the solvent. 

The motor protein jumps at times tj from Uj to nj±d with transition rates k^{nj^ x), 
hydrolyzing (or synthesizing) one ATP molecule per jump which corresponds to tight 
mechanochemical coupling. In this minimal model, we take into account only one 
chemical state. The transition rates have to fulfill a local detailed balance condition 
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of the form 



^ _ , , = exp[(A/i - V{n + d,x) + V{n, x))/k^T] (3) 
k [n + X) 

where we assume that the jumps of the motor protein take place instantaneously. The 
free energy change of the solvent 

A/x = /xatp - Madp - MPi (4) 

is associated with ATP turnover. Implementing mass action law kinetics and the 
concept of a barrier in the potential of mean force for the unresolved chemical steps, the 
individual rates become 

k^{n,x) = k^^exp[AfiATp/kBT]exp[-^ J -^^^^ dn] 

n 

= fc"^exp[A/iATp/fcBr] exp[{-i^d^el/2 - i^{n - x)de+)/k^T] (5) 

and 

n-dO- 

fc-(n,x) = fc^<iexp[(A/iADP + A^pJ/A:Br]exp[-^ j <\n] 



n 



= fc'^exp[(A/XADP + A/xp.)/A;bT] exp[{-Kd^9^/2 + K{n - x)de_)/k^T] (6) 

with 

A/i^ = jii- jf^ = kBT\ii{ci/c^) (7) 

and Ci the concentrations of the nucleotides {i = ATP, ADP, Pi). Here, the transition 
rate k^^ applies to equilibrium concentrations of nucleotides. The load sharing factors 9^ 
and with 9^ + 9_ = 1, depend on the unresolved shape of the free-energy landscape 
of the motor protein [45] . 

2.2. Fokker-Planck equation 

The transition rates, as well as the force the motor protein exerts on the probe, depend 
on the distance 

y{t) = n{t) - x{t) (8) 

between motor and probe. The corresponding probability density p{y^ t) obeys the 
Fokker-Planck-type equation 

dtp{y, t) = k^{y- d) p{y - d,t) + k'iy + d) p{y + d,t) - {k^{y) + k-{y))p{y, t) 

+ dyii^y - fe.)p{y,t) + kBTdyp{y,t))/^, (9) 

which contains in the first line the contributions from the transitions of the motor protein 
and in the second line drift and diffusion of the probe particle. 

For large t and constant nucleotide concentrations, the system reaches a stationary 
state with time independent p^{y) and constant mean velocity (n) = (x) = v with 

(X) = («(?/) - /ex)/7 (10) 
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and 

{h)=d{k+{y)-k-{y)) (11) 

where (...) denotes the average using the stationary distribution p^{y)^ which, however, 
cannot be determined analyticaUy. 

3. Efficiencies 

3.1. First law: single trajectory 

FoUowing the concept of stochastic thermodynamics [46l [47] one can assign a first law 
on the level of a single trajectory. If the probe moves a small distance Ax, the first law 
becomes 

dV 

Agp = (-— - U)Ax = {^y - /ex)Ax (12) 

where Agp is the heat dissipated by the probe, fe^Ax is the work against the external 
force and {dxV)Ax the change of the internal energy of the spring due to the motion of 
only the probe. A jump of the motor protein gives rise to a first law in the form of |^48j 

= AV + A£;soi + A^M 
= AV-Afi + TASsoi + A^M (13) 

without a contribution of the internal energy of the motor as its internal energy does 
not change in the one-state model. The change of the internal energy of the spring is 
given by 

AV = V{n ±d,x)- V{n, x) (14) 

where the sign depends on the the direction of the jump. Due to ATP turnover, the 
internal energy of the solution changes by AE'soi = — A/x + TAaSsoI, where AS'soi is the 
change of the entropy of the solution. The heat dissipated by the motor protein in this 
transition is denoted by A^m- 

3.2. First law: ensemble average 

On average, the chemical energy gained from ATP consumption that involves changes 
of the entropy of the solvent will be dissipated as heat Q in the environment and/or is 
delivered as work against the external force. In the stationary state, the internal energy 
of the spring is constant on average. Taking the average rates of f fT2l ) and f fTSl ) and 
summing the two contributions, this first-law condition can be expressed as 

A/x = (5p + (5m + TSsoi + fe^v (15) 
where the dot denotes a rate and 

A/x = -(Fsoi) = A/x^/rf (16) 
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is the rate of free energy consumption. The rate of dissipated heat Q = Qp + Qm has 
two contributions. First, the heat flow through the motor protein is given by 

Qm = (<?m) = Afi-Vn- TSsoh (17) 

representing the fact that while jumping, the motor protein uses free energy from the 
hydrolysis to load the spring which corresponds to a change of the internal energy of 
the spring 14 with 

/oo 
p'{y)[k^iy)iViy + d)- V{y)) + - d) - Viy))] dy (18) 

-OO 

= ^{k^iy) + k-iy)) + Kd{y{k+{y) - k-{y))). (19) 

The energy thus stored in the spring is then dissipated by the probe whose heat flow is 
given by 

Qp = (gp) = {{Ky - UHy)) (20) 

where 

u{y) = {{Ky - U) + k^Tdy lnf{y))/r (21) 

is the local mean velocity of the probe for a given y ^91 [50] which corresponds to the 
current arising from the motion of only the probe in ([9j). 

3.3. Three different efficiencies 

We will now focus on three diflFerent deflnitions of efliciency that have been proposed 
for motor proteins. 

In the absence of an external force (/ex = 0), one can compare the energy that the 
motor protein transfers to the spring, with its available chemical energy A/x. From 
( ITSl ) and f fTTI ) it follows that Vn = Op- The ratio of on average dissipated heat through 
the probe and available free energy 

Qf _ dK{yv) 

was proposed as deflnition of efliciency [33]. We will see below that t^q is not bounded 
by 1, as it has been anticipated earlier [161 ElJ? ^^id therefore we will call it a pseudo 
efliciency. A second type of efficiency is the Stokes efliciency, 

that compares the mean drag force the probe feels with the available chemical force. 
In contrast to t^q, 77s is bounded by 1 [16j. If the motor protein exerted a constant force 
on the probe, the Stokes efliciency would be equal to the pseudo efl&ciency 77Q because 
in this case the average heat dissipated by the probe is the mean drag force times d. 
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Finally, in the presence of an external force acting on the probe, the thermodynamic 
efficiency of the system is the ratio between mechanical work delivered to the external 
force and available free energy [llOj 

For /ex 7^ 0, the pseudo efficiency t^q can be defined as 

= • (25) 

4. Gaussian approximation 

4.1. Derivation 

For a comparison with the simulations and in order to gain more analytical insights, it 
will be convenient to have a simple approximation for the stationary distribution p^{y). 
For a Gaussian probability distribution 

the free parameters y for the mean and for the variance can be determined by 
requiring that the time-derivative of these quantities as calculated with the Fokker- 
Planck equation ([9j) vanishes in the steady state. These conditions result in the following 
two equations for y and 

{Ky-U)/^ = d{k+-~k-) (27) 

and 

(/^a^ + i^y^ - Uy - ksT)/^ = d[{y - kdO+a'^ /k^T)k^ - {y + kdO^a'^ /k^T)k-] 

+ d\k^ + k-)/2, (28) 
where we have introduced the average jump rates 

/oo 
-00 

= A;"^exp[A^ATp/A:Br - Kd^dlik^T - Ka'^)/2{kBTf - ndd+y/k^T] (29) 
k = k {y)p {y)dy 

J —OO 

= fc"^exp[(A/XADP + A/xfJ/^bT - t^d^elik^T - Ka'^) / 2{kBTf + i^de_y/k^T]. (30) 
These equations can easily be solved numerically. 

4>2. Limits A/x and A/x oo 

Close to chemical equilibrium, i.e.. A// = 0, and for /ex = 0, we expand y and /^a^ — /cbT 
up to first order in A// and find 

y ^ AAfi + A{9^ - e^fAjji (31) 
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and 

^a^ - k^T ^ B{e+ - e_)Afi. (32) 

The coefficients A and B obtained by solving the first order of f l27l ) and f l28l ) are too 
long to be shown here. 

In the limit A/x oc and /ex = 0, we obtain for y and /^a^ — /cbT 

y^^- CAfi exp[- A/x^./^bT] (33) 

and 

i^a^ - k^T ^ DAfi exp[- A/x^./^bT] (34) 

as long as ^_ > 0. The coefficients 

C = 1 + ^d'{0+ - 0-?/^ .35. 

7K#fceq exp[(A^ADP + A^pJ/^bT] 

and 

^ ~27rf2fceqexp[(A/XADP + A/ipJ/^BT] ^^^^ 

are obtained by solving ( l27l l and ( l28l l to first and second order in A/x. 
4^3. Efficiencies 

Within this Gaussian approximation, the average heat fiow through the probe as given 
by ( !20l) is calculated using the local mean velocity f l2T]) 

uiy) = (Ky - U)h - kBTiv - y)/{ja'). (37) 
The average over y can now be performed leading to 

Qp = (K^a' + - UbTk - 2KUy + f^)/^. (38) 
This expression is used to determine rjq in this approximation as 

na^ - k^T + - f^^y 
* = * A,(.y-U <3»> 

with y and being the solution of ( 1271) and (l28l) for given A/x and fc^^. 
For small A/x, using f lHTI) and f l32l) . ryg takes the form 

If ^+ 7^ T^Q diverges for vanishing A/i. For 9^ > tjq can become negative due to 
those jumps of the motor protein that occur when the previous diffusion of the probe 
has resulted in y < — 0.5rf. Then, the energy stored in the spring is dissipated by the 
motor protein during jumping. 

In the limit of large A/x, we use f l34l ) and f l35l ) to obtain 
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Table 1. Values of the model parameters used for the simulation and the Gaussian 
approximation. 

0.407 40 3 • 10^ 



which approaches 1. 

The Stokes efficiency in the Gaussian approximation without an external force is 
simply given by 

For /^cr^ > /cbT, which is the case for 9^ < 0.5, the Stokes efficiency is always smaller 
than rjQ. For vanishing A/x, rjs approaches a finite value, rjs ^ dK.A + dK.A{9^ — ^_)^, 
while for A/i ^ oc it also converges to 1. 



5. Results 



In this section, we study the three efficiencies for our hybrid model as functions of the 
chemical energy A/i, the absolute concentrations of the nucleodides, i.e. fc^^, the external 
force /ex and the load sharing factor 9^. The data are obtained from simulations, using 
a GiUespie algorithm [52j similar to [37j with the motion of the probe being spatially 
discretized in steps of Ax = d/lOOO, and compared with the Gaussian approximation. 

We use model parameters as given in table [H which are motivated by experimental 
results for the Fi^ATPase as described in section [6] below. The load sharing factor 9^ 
remains as a free parameter. 

Simulated trajectories with the same nucleotide concentrations as used in the 
experiment [33j are shown in figure [2l In the presence of low nucleotide concentrations 
only few backward jumps of the motor protein take place and the trajectory of the 
probe shows an almost staircase like form. For high nucleotide concentrations, following 
a forward step the motor often performs a backward jump. Such a sequence of two 
jumps is not necessarily visible in the trajectory of the probe which remains almost 
linear. 



5.1. Pseudo efficiency rjq 

We will first investigate the pseudo efficiency rjq as a function of A/x, fc^^ and We 
extract (jp from the numerical data by averaging over one sufficiently long trajectory. 
The results are shown in figure [3l The most striking fact of these data is the observation 
that T^Q is larger than 1 for small enough A/x and 9j^ which shows up in the Gaussian 
approximation as well. This effect can be understood as follows. In a jump, the motor 
protein can take heat from the solution in order to change the internal energy of the 
spring by an amount larger than A/i. If, subsequently, the probe dissipates this internal 
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Figure 2. Trajectories of the molecular motor (black) and the attached probe particle 
(green) obtained from the simualtion (without external force). In the presence of low 
nucleotide concentrations (I) the trajectory of the probe exhibits a more stepwise form 
while it becomes almost linear for high nucleotide concentrations (III). The parameters 
are (9+ = 0.1, (I): catp = 0.4 /iM, cadp = 0.4 /iM, cp. = 1 mM, (II): catp = 2/iM, 
CADP = 2/iM, Cp. = 1 mM, (III): catp = 100 /iM, cadp = 100 /iM, cp. = 1 mM. For ah 
parameter sets we have A/i = 19.14 /cbT. The nucleotide concentrations are the same 
as used in the experiment [33] shown in figure [8] with the same labelling (I-III) below. 

energy of the spring as heat back into the environment, (jp can indeed become larger 
than A/i without any violation of the second law. Using the obtained parameter for the 
spring constant the motor protein transfers 20 /cbT to the spring if it starts the jump 
from the minimum of the harmonic potential. For small values of the forward jump 
rate of the motor protein depends only weakly on the current position of the probe as 
shown in figure [H Therefore, jumps will occur even if the associated change of internal 
energy of the spring, AV^ is larger than A/x. For rather small k^^^ backward jumps 
are rare and the probe relaxes to the potential minimum between successive forward 
jumps (see data set I in figure [2]). This leads to 14 > on average and hence to 
r]Q > 1 for A/i considerably smaller than 20 /cbT as shown in figure [3l As the value of 9^ 
increases, tjq decreases because the forward jumps of the motor protein are suppressed. 
On average, in this case the motor protein jumps only if the probe has diffused forward 
and exerts a pulling force on the motor through the spring. 

Increasing the absolute concentrations of the nucleotides, i.e., increasing fc^^, results 
in more forward but also more backward jumps, which can be seen for data set III in 
figure [21 For small A/x, the occasional backward jumps follow especially those forward 
jumps for which the change of internal energy of the spring has been larger than A/x, 
leading to a smaller tjq. 

In the limit of large A/i, the motor protein jumps even when the spring is previously 
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Figure 3. Pseudo efficiency t^q from tiie simulation, (a) and (b); and witiiin tiie 
Gaussian approximation, (c) and (d). (a) and (c): 77Q as a function of Aja for different 
values of the load sharing factor 0-^ and fixed k^^ = 10~^s~^. (b) and (d): t^q as 
a function of Aja for different values of k^^ with fixed ^+ = 0.1. The remaining 
parameters are hi = A{)k^T/dP^ 7 = ^A{)lk^T^/dP. In the simulation, the error is of 
the order of the symbol size. 



stretched which can result in changes of the internal energy of the spring by an amount 
larger than 20 k^T. The coupling between the motor protein and the probe induces 
a balancing effect between the forward motion of the motor protein and the drag of 
the probe maintaining a typical V that turns out to be approximately A/i, leading to 
- 1- 



5.2. Stokes efficiency r]s 

We also obtain the Stokes efficiency f l23l) from the simulated trajectories and the 
Gaussian approximation as shown in figure El Characteristically, starting close to for 
small A/x, Tys monotonically increases with A/i reaching 1 for A/x ^ oc. For small A//, 
the trajectory of the probe shows a staircase form with small average velocity leading 
to small values of the Stokes efficiency in contrast to values of the pseudo efficiency 
rjQ > 1. For large A/x, the probe does not relax to the potential minimum between 
consecutive jumps resulting in a more linear trajectory of the probe as if it was exposed 
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Figure 4. Probability distribution p(^)|jump of the distance y just before a jump of the 
motor protein for several values of 6>+ at A/i = 13/cbT and k^^ = 10~^s~^. For small 
the forward jumps of the motor protein are almost independent of the position of 
the probe resulting in a peak at whereas for larger ^+ the peak clearly shifts to 
y < implying that the motor protein prefers to jump when the probe has diffused 
ahead. The peaks around y = 1 indicate backward jumps which take place more often 
in the case of small ^+ when the backward rate is more sensitive to the position of the 
probe. 

to an almost constant force. In this limit of an almost linear motion of the probe, the 
pseudo efficiency becomes the Stokes efficiency. As rjs is bounded by 1, ryg can not reach 
values larger than 1 in this limit either. 

Increasing the load sharing factor 9^ results in decreasing average velocities. 
Therefore, the Stokes efficiency also decreases which can be seen in figures [5] (a) and 
(c). With increasing absolute concentrations of nucleotides, i.e., with increasing /c^^, the 
average velocity and therefore also the Stokes efficiency at fixed A/x increases as shown 
in figures [5] (b) and (d). 

5.3. Thermodynamic efficiency tjt 

The thermodynamic efficiency of the system can be studied only if an external force is 
applied to the probe. As an illustrative example, for fixed Afi = IO/cbT, we examine 
the thermodynamic efficiency in the presence of external forces smaller than the stall 
force as shown in figure [6l The thermodynamic efficiency increases linearly with /ex 
and reaches 1 at the stall force which is possible only due to the tight mechanochemical 
coupling in this model. In figure [6] we also show the pseudo efficiency t^q as defined in 
( l25l) in the presence of external forces. While t^q is almost independent of the external 
force, the contribution from the dissipated heat, (jp/A/x, decreases linearly with /ex and 
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Figure 5. (a) Stokes efficiency r]s from the simulation, (a) and (b); and within the 
Gaussian approximation, (c) and (d). The data is obtained from the same trajectories 
used to obtain 77Q in figure [3] (a) and (b). (a) and (c): rjs as a function of Aja for 
different values of the load sharing factor 0-^ and fixed k^^ = 10~^s~^. (b) and (d): 77s 
as a function of A/i for different values of k^^ for fixed ^+ = 0.1. 



reaches zero at the stall force. 

At stall conditions, the work corresponding to the stall force refers to the maximum 
work the motor protein can convert on average. In the simulation, we find that applying 
/ex = Afi/d generates a diflFusive motion of the motor protein with v 0^ and Qp 
for various values of 9^ and A/i, including the ones with t^q > 1 and t^q < 1. This 
observation implies that the motor protein seems in principle to be able to convert the 
full Afi into extractable work. In our model, where the motor interacts with the external 
force only via the spring, this result is not trivial, as it would have been if we had applied 
/ex directly to the motor. Under these conditions, p^{y) is Gaussian with {y) = A/i/ (dk) 
and the same variance as the Boltzmann-distribution in equilibrium, = kpT/n. 

Within the Gaussian approximation we can insert /ex = A/x/rf in f l27l) and f l28l) . 
With /^cr^ = /cbT, y = A/i/((i/^), i.e., = is a solution for /ex = A/x/rf, implying that 
also in the Gaussian approximation the motor protein is able to convert the full A/x into 
extractable work for any values of the load sharing factors given that 9^ + 9^ = 1. 
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Figure 6. Pseudo efficiency t^q (red and black dots), thermodynamic efficiency r]T 
(green and blue squares) and dissipated heat AQp per A/i (cyan and purple diamonds) 
as functions of /ex for fixed A/i = 19teT, 6>+ = 0.1 and fc^^ = IQ-^s-^ The plot 
contains data from the simulation (S) as well as from the Gaussian approximation 
(G). 



6. Case study: Fi-ATPase 

In this section, we apply our hybrid model to the Fi-ATPase and compare the 
simulations with recent experimental data f33l 1341 [35] . 

6.1. Model parameters 

For a quantitative comparison we have to map the rotary motion of the Fi-ATPase to 
our linear model and determine the model parameters. In our model, one jump of the 
motor protein covering a distance d corresponds to a rotation of the 7 shaft of 120°. 
Using large probe particles like polystyrene beads or act in filaments, the substeps in one 
120° rotation are not resolved experimentally. Therefore, we will omit the substeps here, 
too. We assume that the temperature of the solution is T 24°C and that the probe 
consists of two beads of diameter 287nm [33]. The friction coefficient of the probe can 
be calculated using the formula for the rotational frictional coefficient F from pBl [37] 
with the viscosity of water {rj O.OOlNs/m^). The frictional torque N = Tip acting on 
the probe with angular velocity (f corresponds to a frictional force 
F 

/fr = — X = 7X (43) 

acting on the probe at distance r from the 7 shaft. Within one 120° rotation, the probe 
at distance r covers d = 27rr/3. For the linear model, the friction coefficient 7 can be 
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Figure 7. Comparison of the mean velocities observed experimentally [33 (red 
squares) and in the simulation for several load sharing factors ^+ (black dots, green 
diamonds and blue triangles). The labelling I, II, III, V refers to the corresponding 
parameter sets in figure [H 



calculated as 

^ = ;j = iwr 

leading to 7 = 0.407 fcBTs/rf^. 

Following the mass action law assumption, the equilibrium transition rate fc^^ is 
supposed to depend linearly on the concentrations of nucleotides in the solvent. For 
low ATP concentrations (catp — 10~^M), the mean velocity of the motor protein is 
dominated by the rate of ATP binding. In the one-step model this feature holds for 
all concentrations. Therefore we choose k^^ to be the experimentally determined rate 
of ATP binding fc^^ 3 • 10^M~^s~^c^^p[42j. For known nonequilibrium concentrations 
of nucleotides like in the experiments, the structure of the transition rates ([5j) and ([6]) 
leaves the choice of the equilibrium concentrations arbitrary as long as they obey 

eq 1 

^ 4.89 ■ lO-'i- (45) 

eq eq ly/r V / 

^ADP^Pi 

for T = 23°C and pH 7 [37j. For given k^^ and A/x, one possible choice of the 
nonequilibrium concentrations of nucleotides is cadp = c^dp' ^Pi ^ ^atp = 

c^p exp[A/x//cBr] which was used for the simulation and the Gaussian approximation. 

In order to determine the spring constant and the load sharing factor 9^ we use 
both the experimental data of the mean velocities [33j and the histogram of the angular 
position of the probe at a jump [35j. While both data sets depend on both parameters, 
the velocity, especially for large /c^^, is more sensitive to whereas the peak position of 
the histogram mainly depends on 9^. Therefore, we primarily use the velocity data to 
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Figure 8. Average heat Qp released through the probe (green and cyan) and work 
W = fexd against the external force (blue) compared to the available free energy per 
step A/i (red line). The dissipated heat of the probe is split into two contributions 
Qs and Qy according to the two terms of the Harada-Sasa relation ()A.3p . The 
contribution from the linear motion with constant mean velocity, Qs (cyan), appears 
in the numerator of the Stokes efficiency while Qv (green) is the contribution due to 
the non-uniform jumping motion of the motor protein. In each of the five parameter 
sets labelled by I-V, the left and the central bar represent results from the simulation 
for 0-^ = 0.1 and ^+ = 0.01, respectively, while the right bar shows the experimental 
results and error bars from [33]. The following parameters were used in the five cases: 



(I) CATP = 0.4 /iM, CADP 

A/i = 19.14 ^bT; (II) CATP 



and 



= 0.4 /iM, cp. = 1 mM, i.e, /c^^ = 5.87 • 10" 
2/iM, CADP = 2/iM, Cp. = 1 mM, i.e, fe^^ = 2.93-10-'^ s" 



and A/i = 19.14 ^bT; (III) catp = 100 /iM, cadp = 100 /iM, cp. = 1 mM, i.e, 
A;^^ = 1.47- 10-^s-i and A/i = 19.14 /cbT; (IV) catp = 2/iM, cadp = 2/iM, 
Cp. = 1 mM, /ex = 9.27 ksT/d, i.e, fc^^ = 2.93 • 10"^ s"^ and A/i = 19.14 ^bT; 
(V) CATP = 2/iM, CADP = 0.5 /iM, Cp. = 0.5 /iM, i.e, A:^^ = 3.67- 10"^^ s"^ and 
A/i = 28.12 A:bT; 



fit and determine the load sharing factor 9^ by comparing the peak position of the 
experimental histogram [35j with the left peak position of the corresponding histograms 
obtained by our simulation as the ones shown in figure [H 

As a result, we obtain = 40 ± bk^T/d^ and a value of 0+ in the range 
0.1 < ^+ < 0.3. In figure [Zl we show how for this value of n changing the load sharing 
factor affects the mean velocity for which we get the best overall agreement for 0+ = 0.1. 
For later purposes, we also include data for 0+ = 0.01. 

6.2. Comparison of efficiencies with experimental data 

Experimentally, the heat flow of the probe is determined using the Harada-Sasa relation 
[53] . In the appendix, we show that this heat flow is equal to (jp as defined in ( l20l) . 
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In figure [HI we plot the average heat released through the probe per step, Qp, plus the 
work against the external force, obtained by the simulation for = AQk^T/d^ 
and = 0.1 and compare it with the experimental results [33j. We find quite 
good agreement between theory and experiment for the parameter sets I-IV where 
either the maximum deviation is 15% (II- IV) or our theoretical value is included in 
the experimental error range (I). As an aside, we note that for the parameter sets 
I-III (without external force) also the simulated mean velocities coincide well with 
the experimental values with a maximum deviation of 10% as shown in figure [3 For 
illustrative purposes, we also plot Qp plus W for 9^ = 0.01 which shows better agreement 
with the experimental data (but is not consistent with the range of 9^ obtained in section 

EH- 

Discrepancies between our theory and the experiment are visible in figures [71 and [8] 
where for parameter set V both the average velocity and the pseudo efficiency deviate 
significantly from the experimental values for = AQk^T/d^ and 0+ = 0.1. For 
A/x = 28.12/cBr corresponding to the data set V in figure [HI the probe just reaches 
the potential minimum between consecutive jumps of the motor protein. Therefore on 
average at most 20/cbT can be transferred to the spring leading to t^q c::^ 0.7, which is 
less than the experimental value. This discrepancy is most likely due to the fact that 
we have omitted substeps in our model. In the simulation, the average velocity then 
does not show saturation as one would expect it to result from the hydrolysis step [42j 
which should be experimentally observable at the higher concentrations used in [33j. 

The confinement of ^+ to the range 0.1 < 6^+ < 0.3 implies on the one hand that 
the potential of mean force of the motor protein should be asymmetric and on the other 
hand that asymmetric potentials with a barrier state close to the initial state seem to 
enhance the ability of the motor protein to perform work on the spring, in accordance 
with [54]. If 9j^ was larger, t^q would decrease and the experimentally determined values 
of T^Q would not be reached in the simulation. If ^+ was smaller, t^q would approach the 
experimental values better, however, the distribution of the position of the probe just 
before a jump as shown in figure [4] would then no longer coincide with the experimentally 
observed distribution (see [35j). 

Information about the thermodynamic efficiency of the motor protein can be gained 
by applying an external force at the probe. In our simulation, the stall force is found 
to be /ex = A/x/rf, implying that the motor protein is able to convert the full A/x into 
extractable work without dissipation in accordance with the experiments performed in 

7. Conclusion 

In summary, we have discussed a simple generic model which includes the elastic linker 
between the probe particle and the molecular motor. Properties of the motor become 
typically accessible only through the observation of the motion of the probe. We have 
then focussed on discussing three types of efficencies within this model using both 



Efficiencies of a molecular motor: A generic hybrid model applied to the Fi^ATPase 18 



simulations and a Gaussian approximation to the stationary distribution for the distance 
between motor and probe. The genuine thermodynamic efficiency is non-zero only if 
an external force is applied to the probe. The Stokes efficiency deviates from 1 due 
to the discrete nature of the motor steps which become less relevant with increasing 
ATP concentration. A pseudo efficiency measuring how much of the free energy of ATP 
hydrolysis ends up in loading the elastic element can even become larger than 1 close 
to equilibrium and for a barrier state close to the initial state. 

Applying this minimal model to recent experimental data for the Fi-ATPase we find 
overall good agreement except for those parameters where especially the Pi concentration 
is very small. In general, one should consider ATP binding and ATP hydrolysis as two 
separate steps. Such a refinement as well as a splitting of the 120° rotation into two steps 
of 90° and 30° as experimentally observed using much smaller probe particles does not 
pose new conceptual challenges to the present framework and will be pursued elsewhere. 

Appendix A. Equivalence of heat flow (jp with the one inferred from the 
Harada— Sasa relation 

Experimentally, the heat fiow caused by the probe has been inferred from measuring the 
autocorrelation function Cx{t) = {x{t + T)x{t)) — v'^ and the linear response function 

^x(r)_ ^^^^^ (A.l) 

of the velocity of the probe to a small external perturbation h{t) of the probe within 
the steady state [33j . The heat fiow is then given by an equality derived by Harada and 
Sasa ^53j 

gHS = ^^2^^ r ^[c. _ 2/cBTRe (i?,(c.))] 

= W + 7[Cx(0) - 2kBTRi{0)] (A.2) 

= + Qv, (A.3) 

with Cx and Rx being the Fourier transforms of Cx and Rx. Using a path weight 
approach described in [55j applied to our system, the response function follows as 

R^ir) = ^{^it + ^Mt) - i^yit) - /ex)/7])- (A.4) 
Inserting Cx and this Rx into f lA.2[ ). one immediately finds 

= _ /^j) (A.5) 

which is equal to (jp in f l20l ). 
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